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Topological susceptibility in Yang-Mills theory 
in the vacuum correlator method 
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We calculate the topological susceptibility of the Yang-Mills vacuum using the field correlator method. 
Our estimate for the SU(3) gauge group, y}^ = 196(7) MeV, is in a very good agreement with the results of 
recent numerical simulations of the Yang-Mills theory on the lattice. 
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The topology of the QCD vacuum and the hadron 
phenomenology are intimately related. The celebrated 
example is given by the Witten-Veneziano formula pQ 
which provides a solution to the "U(l) problem" by uti- 
lizing the explicit breaking of the U(l) axial symmetry 
in the large- N QCD. In a simplest case of Nf massless 
quarks the formula relates the mass of the rf meson, 
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to the non-perturbative susceptibility 
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of the topological charge 
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calculated in the pure SU(N) Yang- Mills theory and 
to the pion decay constant (note that (Q) = due 
to CP-invariance of the Yang-Mills vacuum) . An exact 
value of the susceptibility x is the subject of intensive 
numerical studies in the lattice gauge theories [H [3] . 
The local density q(x) of the topological charge is 
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where F^ — F^ v T a is the non-Abelian field strength 
tensor with the components 
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gf 



abc a b 
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and T a , a — 1, . . . , N 2 — 1 are the generators of the 
SU (N) gauge group normalized by the standard rela- 
tion, Tr T a T b = 5 ab /2, and f abc = —2i Tr[T a , T b ]T c is 
the totally asymmetric structure constant of the group. 



In this paper we estimate the topological suscepti- 
bility ([2]) in the formalism of gauge-invariant non-local 
correlators _T which is reviewed in Ref. [5J [5] . We work 
in Euclidean space-time suitable for comparison of our 
results with the ones obtained in numerical simulations 
of the lattice Yang-Mills theory. A basic element of the 
field correlator method is a non-local quantity, 



G MV (x', x ) = ${x ; x)F MV (x)$(x', x ) , 



(6) 



where $>(x;xo) = <f>'(xQ,x) is an ordered exponent 
of the integral over the gauge field = T a A^ (the 
Schwinger line). The integration is taken along the ori- 
ented path C X-XQ which stretches between the points Xq 
and x of the space-time, 



&(x; y) = Vexp 
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A M (z) dz^ 
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where V denotes path ordering. Under the SU(N) 
gauge transformation fl the quantity ([7]) changes as 

*(x;v)^n(x)$(x;y)tf(y). 

Due to this property the non-local object ^ transforms 
essentially locally, 

G fiv (x, x ) -> rt(x )G fll ,(x, x )& (x ) , 

making it possible to construct various (non-local) 
gauge invariant objects of the general form 

D$...v n 00 "^[^Mifi^i^o) ■■■G tlnVn (x n ;x )] . (8) 

One may think about this quantity as of the n-point cor- 
relation function of the field strength tensors © covari- 
antly shifted from the "reference point" xo to the points 
Xi, i = 1, . . . ,n along appropriate paths. The simplest 
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nontrivial quantity of the kind (jHJ) is the two-point cor- 
relator (hereafter we omit the superscript n = 2) 

g 2 

D^ a /3(xi,x 2 ;x ) = —Tr[Gp V (xi,xo)Gap(z2',%o)\ -(9) 

One may choose the paths connecting the points Xi, 
i = 0, 1, 2 as straight lines and put the point :ro exactly 
on a straight line between the points x\ and xi. Then 
the correlator (J9j) becomes a function of a single variable 
z = x\ — x 2 . 

At zero temperature the correlator (03) can be pa- 
rameterized by two scalar functions D and D\ 
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(z) = (S^fi - 5tf5 va ) D{z 2 ) (10) 
zp8 va ) -{n<-+v) Di(z 2 ) . 
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Perturbatively, both functions are divergent as 1/|^| 4 at 
short distances, \z\ — > 0. However, besides the pertur- 
bative part the structure functions possess also a non- 
perturbative part which is relevant to the long-distance 
physics. The lattice simulations [Jj [8] indicate that the 
structure functions can be well described as follows 



D(z 2 ) = A e-^/ T * 
Di(z 2 ) = A ie ~^/ T * 



e -\z\/X 
P-I^l/A 
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where the first terms in both expressions correspond to a 
non-perturbative contribution while the last two terms 
contain the perturbative parts. The functions D and 
D\ are parameterized by the correlation lengths T g and 
A, and by the prefactors Aj and b{, i = 0, 1. 

On general grounds one may expect that the corre- 
lators ([8]) should depend not only on the points Xi, but 
also on the auxiliary variables such as xq and shapes 
of the paths C Xn , X(l entering the phase factors ([7]) used 
in Eq. ([6]). Perturbative corrections coming from the 
Schwinger lines may provide an additional contribution 
to the correlation lengths. However, a physical non- 
perturbative observable - if expressed in terms of the 
correlators (J5J only - should be independent of both the 
shape of the Schwinger lines and the ultraviolet cutoff. 

In what follows we consider the quadratic correla- 
tor ([9]) for which various numerical estimations are avail- 
able [7J [51 01] • The numerical results were obtained in 
lattice simulations with the use of the so-called cooling- 
method which makes it possible to get rid of short-range 
lattice artifacts while leaving the long-range physics in- 
tact. The cooling removes ultraviolet suppression of 
long Schwinger lines <&(x,y), which may originate from 
perturbative fluctuations. As a result of the (soft) cool- 
ing, the correlation lengths T g and A in the correla- 
tors (fTTj) do not depend on the ultraviolet cutoff both in 



pure Yang-Mills theory [7J and in QCD with dynamical 
quarks [7] . Similar results were obtained in Ref. 03 with 
the use of the so called renormalization group smoothing 
which is an alternative to the cooling. 

Another subtlety of the lattice results is the depen- 
dence of the correlators on the shape of the Schwinger 
line connecting the points xq, x\ and x 2 in the defi- 
nition of the correlator (03). It turns out that the de- 
pendence of the prefactors Ai and i = 0, 1 is rather 
strong, while the correlation lengths T g and A are insen- 
sitive with respect to the variations of the shape of the 
line [10]. Thus our results should be understood with 
the prescription that the perturbative ultraviolet correc- 
tions are subtracted and the Schwinger line itself is cho- 
sen to take a most natural shape forming straight path 
stretched between the points x\ and x 2 in the definition 
of the correlation function (03). In this prescription the 
correlator is clearly independent of the reference point 
xq provided it is located on the Schwinger line stretched 
between the points X\ and x 2 . 

It is argued [5; that in the Yang-Mills theory and in 
QCD the dominant contribution to various observables 
is given by the lowest bilocal correlator. This obser- 
vation is often call as the "Gaussian dominance". The 
validity of the Gaussian approximation is supported, for 
example, by the argument, that the contribution of the 
bilocal correlations to the interactions between colored 
particles dominates the contribution coming from the 
higher-order correlators. The dominance could be inter- 
preted either in stochastic picture (any n > 2 correlator 
contributes much less then the binomial one) or in the 
coherent picture which postulates that there is a strong 
cancelation between all n > 2 correlator so that the 
result is given only by the binomial contribution [llj . 
We take an advantage of the stochastic scenario assum- 
ing that all higher order correlators are suppressed with 
respect to the leading Gaussian contribution. 

In the Gaussian vacuum all n-point correlators of 
field strengths are assumed to be factorized into the 
bilocal correlators. In particular, this means that all 
odd correlators are zero in the Gaussian approximation. 
The factorization of even correlation functions goes ac- 
cording the scheme [S]: 

(G Ql (l)G Q2 (2) . . . G a2 "(2n)) - {S aia2 . . J" 20 - 1 " 3 -) 
x (Tr G(1)G(2)) . . . (Tr G(2n - l)G(2n)) (12) 
+all permutations , 

where the color arguments of the field strength tensors 
are written in a short form. 

The quadratic correlation function ([2]) of the topo- 
logical densities involves the correlator of the four field 
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strengths operators located in two points of the space 
time. A related four-point correlator of field strengths 
can be calculated in the Gaussian approximation ac- 
cording to the factorization scheme described above: 

(K^)F b lS (Q)F^{z)F%{z)) 

= <F tt yO)F 7 & 5 (0)) (F^(z)F d e (z)) 

+{K^)nM) (F b s (0)F d d (z)) 

HF^)F%(z)) (F'MF^iz)) 



(13) 



N 



5 ab S cd D a/3lS (0)D^ p g(0) 



\ n 2 - 1 
+S ac 5 bd D a ^ v [z)D l5pe {z) 
+5 ad 8 bc D aPp6 (z)D^ u {z) 

In this calculation we identified the field strength 
tensors F^ u - used in the definition of the topological 
charge density ([4J - with the covariantly transformed 
field strength tensors G M „, Eq. ©, used in the def- 
inition of the gauge-invariant bilocal correlators ©. 
We used the observation that the four-point correlation 
function (fT3"|) involves only two distinct points x\ = 
and X2 — z. As a result, the expression ([13")) can be fac- 
torized into two bilocal correlation functions with the 
same Schwinger lines. In the axial gauge, z^A^z) = 0, 
the Schwinger line ([7]) is equal to unity, and therefore 
one arrives to the identification G = F, which is valid 
only in this gauge. However, since both the topological 
density ([J| and the bilocal correlator © are the gauge- 
invariant quantities, our results - obtained with in the 
axial gauge - must be gauge-independent. 

The correlation function of two topological densi- 
ties dU) is calculated with the help of Eqs. dTUJ) and ([Hi]) : 



1 



32tt 4 N 2 
: {2 [£>(«) + £>i(z) 



D(z)+D 1 (z) 
dDx(z 



+ z u 



dz u 



} ■ (14) 



In a physical language the use of the factorization pre- 
scription (|13[) in the pseudoscalar channel implies that 
we have associated (a tower of) pseudoscalar glueballs - 
which mediate the interaction between the pseudoscalar 
g-probes - by two non-perturbatively dressed gluons. 
For a more appropriate treatment of the glueballs in 
the field correlator formalism see Ref. [12] . 

It is important to notice that the two-point correla- 
tion function (| 14[) in Euclidean space-time must always 
be negative for non-zero z, 



(q(0)q{z))<0 for \z\ > . 



(15) 



due to a reflection positivity property and a pseu- 
doscalar nature of the topological charge [13] . Using the 



parametrization for the bilocal correlators (fTTj) as well as 
the fitting results obtained with the help of the lattice 
simulations [8] [14] one can immediately check that 
the negativity requirement is indeed satisfied for large 
enough distances (typically, for \z\ > 1.5 fm). However, 
at smaller distances the correlator ([14]) becomes positive 
due to large positive contributions coming from the per- 
turbative parts of the functions D(z) and Di(z). These 
parts are given in Eq. (jTTJ) by the terms proportional to 
the factors bo and b±. 

It is important to stress that the positivity of the 
correlator ([14]) (calculated with the parameters taken 
from the lattice data [7] [8] [14] ) has an apparent incon- 
sistency with the asymptotic freedom of the Yang-Mills 
theory. In fact, at small distances the physics must be 
dominated by the perturbative corrections. For exam- 
ple, in the tree order (here ot s — g 2 /(47r)) one has 

D^(z) = 0, D r°(z) = ^.^ljL, (16) 

7T N |z| 4 

which leads to the known negative- valued result, 

3a 2 N 2 - 1 



4tt 6 



< 0. 



(17) 



One can check that perturbative corrections to D(z) 
and Di(z) structure functions [15] do not spoil this re- 
sult. Note that we do not discuss here non-perturbative 
physics at short distances since the corresponding terms 
must anyway be small compared to the leading pertur- 
bative result. 

Moreover, in order for the correlator (|15[) to be neg- 
ative at small distances, the parameters bo and b\ of the 
parametrization (jTTJ) must satisfy the condition 



bo < h. 



(18) 



However, neither of available results of the lattice simu- 
lations [Z1181Q3] of the pure Yang-Mills theory and QCD 
agrees with the requirement (|18p . 

The observed inconsistency with the perturbative 
part of the lattice results and the reflection positivity 
does not however undermine non-perturbative calcula- 
tions utilizing the field correlator method. A plausible 
explanation of the contradiction may be related to the 
cooling procedure used in the lattice simulations. In 
fact, the cooling may affect the short range part of the 
bilocal correlators, leading to a modification of the co- 
efficients 60 and 61 as compared with their values in the 
uncooled vacuum. In Ref. [7] it is clearly demonstrated 
that the field correlator at a fixed short separation, e.g. 
\z\ « 0.3 fm, in the course of the cooling increases by a 
factor of 3 from its initial negative value to a positive 
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value at a plateau. Moreover, the fits of the numerical 
data for the correlators were done at higher separations 
z, and therefore we do not expect that the fitting results 
represent correctly the short distance physics. As a con- 
sequence, there is no contradiction between the lattice 
results [HOE] and the reflection positivity ([15]) . 

Since the aim of this paper is to evaluate the non- 
perturbative contribution to the susceptibility of the 
topological charge, our results should not be affected by 
uncertainties in determination of the perturbative part 
of the bilocal field correlators. 

Another important remark is that in our calcula- 
tions given below we neglect the contact term Q3] in the 
correlator of the local topological densities (jT4j) . This 
term should appear as a singular 5-like function at zero 
separations, z = 0, in order to reconcile the apparent 
positivity of the susceptibility ([2]) with the negativity 
requirement p5[) . Unfortunately, the contact term can- 
not be accounted in our calculations because we are 
constrained (by results of the lattice simulations) to 
use the particular prescription (jlip of the field corre- 
lators (jTUJ) . Moreover, we expect that in the quoted 
lattice calculations this term cannot be accessible any- 
way because these numerical calculations use the cooling 
method which must destroy (or affect substantially) the 
short distance physics. On the other hand, we know 
that the (soft) cooling does not affect the chiral and 
bulk topological properties of the QCD vacuum [IB]. In 
terms of the correlation function (fT4| the latter obser- 
vation may indicate that in the course of the cooling the 
non-perturbative part of the zero-distance singularity is 
not destroyed literally but it is rather shifted towards 
longer distances. If true, this means that we can ne- 
glect the z — singularity provided we use only the 
non-perturbative z ^ part which is obtained with the 
use of the cooling method. With this justification in 
mind we now continue with explicit calculations. 

We use the data of Ref. [14] where the long-range tail 
- left intact by the cooling procedure - was fitted only 
by the non-perturbative terms of Eqs. (|11[) . The fits 
(with vanishing perturbative part, bo = b% = 0) provide 
us with the following results [14] 

L> NP (0) = A = 3.62(19) A| = 0.212(11) GeV 4 , 
D^ P (Q) = Ax = 1.23(7) A| = 0.072(4) GeV 4 , (19) 

' ' 0.222(4) fm, 



Tg 183(3) A 



0.900(14) GeV 



written in various units for the sake of convenience. 
Here the superscript "NP" stands for "non-perturba- 
tive", and Al — 4.92 MeV denotes a lattice renormal- 
ization scale chosen in Ref. [14] , 



Substituting Eq. (TT4"]) into Eq. © and taking non- 
perturbative parts of the correlators (fTT]) we arrive to 
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N 2 



X 



D NP (Q) D Nt, (0) + D^(0) T 4 . (20) 
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lNP, 



NP/ 



Finally, a direct evaluation of the susceptibility ([211]) 
with the help of the non-perturbative part of the field 
correlators gives us: 



1/4 



Xth. 



196(7) MeV. 



(21) 



Here we used the parameters (Til?]) from Ref. [2] as a 
numerical input. 

It is remarkable that our result (|2ip - based on the 
direct evaluation of the topological susceptibility - co- 
incides (within the small error bars) with the numerical 
value 



Xiauico = 193(9) MeV . 



(22) 



obtained recently in lattice simulations of pure SU(3) 
gauge theory [3]. One should notice, however, that the 
value of the error in our estimation (|21[) comes only 
from the numerical error of the lattice data (fl"9)) . Con- 
sequently, the error in Eq. ([2"T]) does not reflect any 
systematic uncertainty which could be related to our 
assumptions implied in the course of the derivation of 

Eq. dsn}. 

In conclusion, we stress that in our calculation of 
the topological susceptibility we have used the partic- 
ular prescriptions related to (i) the renormalization of 
the ultraviolet divergences (self-energy) associated with 
the Schwinger lines, (ii) the path-dependence of the field 
correlation functions, and (iii) the presence of the diver- 
gent contact term in the correlators of the topological 
charge densities. After a proper treatment of these sub- 
tle issues we arrive to the analytical result for the topo- 
logical susceptibility ([20]) . In the case of SU(3) gauge 
theory we get the numerical value (|21[) which is in a very 
good agreement with the recent lattice estimation 
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